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MULTIPLE 
INTEGRALS 


The first seven sections of this chapter develop the double and triple integral. They 
depend on Sections 11.1 and 11.2 on surfaces and continuous functions, but are 
independent of Chapter [0 on vectors. 

Sections 8 through 10 of this chapter discuss the relationship between 
multiple integrals, line integrals, and surface integrals. Chapters 10 on vectors and 11 
on partial derivatives are prerequisites. 


DOUBLE INTEGRALS 


The double integral is the analogue of the single integral (definite integral) suggested 
by Figure 12.1.1. Figure 12.1.1(a) shows the area A bounded by the interval [a, b] and 
the curve y = f(x), and corresponds to the single integral 


A= [ sora 


Figure 12.1.1(b) shows the volume V bounded by the plane region D and the surface 
z = f(x, y), and corresponds to the double integral 


V= {fro dx dy. 


D 


Our development of double integrals will be similar to our development of 
single integrals in Chapter 4. Before going into detail, we give a brief intuitive preview. 

Instead of closed intervals [u, v] in the line, we deal with closed regions D 
in the plane. A volume function for f (x, y) is a function B, which assigns a real number 
B(D) to each closed region D, and has the following two properties: Addition Property 
and Cylinder Property. 
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I(x) 


(a) Area (b) 


Volume 


Figure 12.1.1 


ADDITION PROPERTY 


If D is divided into two regions D, and D, which meet only on a common 
boundary curve, then 


B(D) = B(D,) + B(D,). 


(Intuitively, the volume over D is the sum of the volumes over D, and D,.) 
This property is illustrated in Figure 12.1.2(a). 


CYLINDER PROPERTY 


Let mand M be the minimum and maximum values of f(x, ) on D and let A be 
the area of D. Then 


mA < B(D) < MA. 


(Intuitively, the volume over D is between the volumes of the cylinders over 
D of height m and M. This corresponds to the Rectangle Property for single 
integrals.) 


This property is illustrated in Figure 12.1.2(b). 


We shall see at the end of this section that the double integral 


Jv (x, y) dx dy 


is the unique volume function for a continuous function f(x, )). The double integral 
will be constructed using double Riemann sums, just as the single integral was 
constructed from single Riemann sums. 

We now begin the construction of the double integral, starting with a careful 
discussion of closed regions in the plane. 


A closed region in the (x, y) plane is a set D of real points (x, }) given by 
inequalities 


a, Sx <a, b(x) S y S b,(x), 
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(a) Addition Property (b) Cylinder Property 
Figure 12.1.2 


where b,(x) and b,(x) are continuous and b,(x) < b,(x) for x in [a,,a,] (Figure 
12.1.3(a)). The boundary of D is the set of points in D which are on the curves 


x= ay; X = ap, y = b,(x), y = b,(x). 
The simplest type of closed region is a closed rectangle 
a, SxS a), by Sy by, 
shown in Figure {2.1.3(b). 


Remark In this course we are restricting our attention to a very simple type of closed 
region, sometimes called a basic closed region. In advanced calculus and 
beyond, a much wider class of closed regions is studied. 


An open region is a set of real points defined by strict inequalities of the 


form 
C, <X <p, d,(x) < y < d,(x). 
y y 

by F--- 

bi---} 
+ - 

x ae x 
(a) Closed region (b) Closed rectangle 


Figure 12.1.3 


713 


714 


12 MULTIPLE INTEGRALS 


We shall usually be working with closed regions. So from now on when we 
use the word region alone we mean closed region. 
To simplify our treatment we shall consider only continuous functions. 


PERMANENT ASSUMPTION FOR CHAPTER 12 


Whenever we refer to a function f(x, y) and a region D, we assume that f(x, y) 
is continuous on some open region containing D. 


If f(x, y) = 0 on D, the double integral is intuitively the volume of the solid 
over D between the surfaces z = 0 and z = f(x,y); Le. the solid consisting of all 
points (x, y, z) where (x, y) is in D and 


0<z<f(x,)). 


If f(x, y) < 0 on D the double integral is intuitively the negative of the volume 
of the solid under D between the surfaces z = {(x, ) and z = 0. Thus volumes above 
the plane z = 0 are counted positively and volumes below z = 0 are counted 
negatively (Figure 12.1.4). 


Figure 12.1.4 


We now define the double Riemann sum and use it to give a precise definition 
of the double integral. We first consider the case where D is a rectangle 


a, Sx Sa), b, Sy S by, 


shown in Figure 12.1.5. 

Let Ax and Ay be positive real numbers. We partition the interval [a,, a2] 
into subintervals of length Ax and [b, , b,] into subintervals of length Ay. The partition 
points are 


Xo =4,;, X, =, + AX, X2 =a, + 2AX,...,x, = a, + Ax, 
You by, yy = db) + Ay yy = dy + 2AY,... dp = 1 + PAY 


where X, <4, SX, + Ax, Vp Ss CxS P,P ay. 
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J 
by +----— 
D 
b1 lew Sl | 
| | 
ay a2 i 


Figure 12.1.5 


If Ax and Ay do not evenly divide a, — a, and b, — b,, there will be little 
pieces left over at the end. We have partitioned the rectangle D into Ax by Ay sub- 
rectangles with partition points 


(X;5 Ys O<k<a, O<Il<p, 


as in Figure 12.1.6. 


Figure 12.1.6 


The double Riemann sum for a rectangle D is the sum 
n p 
Dd L(x, y)Ax Ay = Yo DT FO, yy) Ax Ay. 
D k=01=0 


This is the sum of the volume of the rectangular solids with base Ax Ay and height 


I (Xes Yi): 


As we can see from Figure 12.1.7, 
Dd S(, y) Ax Ay 
D 
approximates the volume of the solid over D between z = 0 and z = f(x, y). 
Now let D be a general region 
a, Sx Sa), b,(x) S y S bax). 
The circumscribed rectangle of D is the rectangle 


a, 5x54), B,SySB,, 


715 


716 12 MULTIPLE INTEGRALS 


Zz Double Riemann Sum 


Figure 12.1.7 


where B, = minimum value of b,(x). 


8B, = maximum value of b,(x). 


It is shown in Figure 12.1.8. 


Figure 12.1.8 


Given positive real numbers Ax and A). we partition the circumscribed 
rectangle of D into Ax by Ay subrectangles with partition points 


(Xn. 3) Osksn, O</l<p. 


DEFINITION 


The double Riemann sum over D is defined as the sum of the volumes of the 
rectangular solids with base Ax Ay and height f(X,..,) corresponding to 
partition points (x,, ¥,) which belong to D. In symbols, 


fie WANApe- > )>., (igen aear, 
dD 


(xg.y7) in D 


Notice that in the double Riemann sum over D, we only use partition points 
(x,.4;) which belong to D (Figure 12.1.9). 
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Figure 12.1.9 Double Riemann Sum 


EXAMPLE 1. Find the double Riemann sum 
D> x?y Ax Ay, 
Dy 


where D, is the square 


and 


The partition of D, is shown in Figure 12.1.10 and the values of xy at the 
partition points are shown in the table. 


2 = =i —2 = 2 =, 
x"y Yo = 0 Vi=5 y2=5 V5.5 y4=35 
Xp = 0 0 0 0) 0 0 

pee i _ we =e 4 
X,=4 0 80 80 80 80 

ees a 8 12 16 
x2 =72 0 80 80 80 80 

L | 

per) 19. 18 27 36 

x3 =4 0 80 80 80 80 


The double Riemann sum is 
Yd x?y Ax Ay 
Dy 
=(142434444484 124 16494 18 +27 + 36}45-4-4 = 0.0875. 


A similar computation with Ax = 75, Ay = 7p gives 
YY x?y Ax Ay = 0.12825. 
Di 
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EXAMPLE 2 Find the double Riemann sum 
DD x7y An Ay, 
D, 


where D, is the region 


and 


The circumscribed rectangle of D, is the unit square. The partition and D, 
are shown in Figure 12.1.11 and the partition points which actually belong to 
D, are circled. The table shows the values of x?) at the partition points which 
belong to D,. It is a part of the table from Example 1. 


eo Noe WE be? HE PSs ess 
Xp = 0 | 
: I 
x, = 4 ay 4 | 
7 j 
X2 = 5 36 | 86 | 
| ! i I 
: | : 
x3 = 3 | | 8 1 %8 
The double Riemann sum is 
YS x?y Ax Ay 
D2 
12 27 36,1 1 86 
= + é aie g + = 0.05375. 


80 * 807 80 ~ 807 80” 80/4°S 80-4-5 
A similar computation with Ax = 75, Ay = 7g gives 
>>, x?y Ax Ay = 0.04881, 
D2 


y v 
} 7 
3 ae Geen | 
3 4 + + 
5 
dD, 
5) aaa ESS GE 
1 
5 ; a 5 | 
0 ws ~+ ‘. 
Ji rs ce er nr 
4 4 4 


= 


1 


Figure 12.1.10 Figure 12.1. 
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Given the function f(x, y) and the region D, the double Riemann sum 


LS (x, y) Ax Ay 
D 


is a real function of Ax and Ay. When we replace Ax and Ay by positive infinitesimals 
dx and dy (Figure 12.1.12), we obtain (by the Function Axiom) the infinite double 


Riemann sum 


LLL y) dx dy. 


The infinite double Riemann sum is in general a hyperreal number. Intuitively, it is 
equal to the sum of the volumes of infinitely many rectangular solids of infinitesimal 
base dx dy and height f(x,, y,). The double integral is defined as the standard part of 
the infinite double Riemann sum. The following lemma, based on our Permanent 
Assumption for Chapter 12, shows that this sum has a standard part. 


Infinite Double Riemann Sum 


Figure 12.1.12 


LEMMA 
For any positive infinitesimals dx and dy, the double Riemann sum 


>> J (x, y) dx dy 


is a finite hyperreal number and thus has a standard part. 


We omit the proof, which is similar to the proof that single Riemann sums 
are finite. We are now ready to define the double integral. 


DEFINITION 


Given positive infinitesimals dx and dy, the double integral of a continuous 
function f(x, y) over D is the standard part of the double Riemann sum: 


{ [res y) dx dy = st [Dy slo sdax ay) 
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Here is a list of properties of the double integral. Each property is analogous 
to a property of the single integral given in Chapter Four and has a similar proof. 
INDEPENDENCE OF dx AND ay 


The value of the double integral |\p f(x, y) dx dy does not depend on dx and dy. 
That ts, if dx, d,x, dy, and d,y are positive infinitesimals then 


[[resnacar= [fronacay 
D 


D 


This theorem shows that the value of the double integral depends only on 
the function f and the region D. From now on we shall usually use the simpler notation 
dA = dx dy for the area of an infinitesimal dx by dy rectangle, and 


| [res y)dA for | [re y) dx dy. 
D D 


ADDITION PROPERTY 


Let D be divided into two regions D, and D, which meet only on a common 
boundary as in Figure 12.113. Then 


{ fos yd = | [ fy yd + | i fle, »)dA. 
D D, D2 


Figure 12.1.13 (a) (b) 


Interpreting the double integral as a volume, the Addition Property says that 
the volume of the solid over D is equal to the sum of the volume over D, and the 
volume over D,, as shown in Figure 12.1.14. 

A continuous function z = f(x,y) always has a minimum and maximum 
value on a closed region D. The proof is similar to the one-variable case. 
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Figure 12.1.14 Addition Property 


CYLINDER PROPERTY 


Let mand M be the minimum and maximum values of f(x, y) on D and let A be 
the area of D. Then 


mAs [free y)dA < MA. 


D 


This corresponds to the Rectangle Property for single integrals. The solid 
with base D and constant height m is called the inscribed cylinder, and the solid with 
base D and height M is called the circumscribed cylinder. The inscribed cylinder and 
the circumscribed cylinder are shown in Figure 12.1.15. Intuitively, the volume of a 
cylinder is equal to the area of the base A times the height. Thus the Cylinder Property 
states that the volume of the solid is between the volumes of the inscribed and circum- 
scribed cylinders. 

Here are two consequences of the Cylinder Property. 


circumscribed 
cylinder 


inscribed 
cylinder 


Figure 12.1.15 Cylinder Property 
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COROLLARY 1 


The area of D is equal to the double integral of the constant function | over D, 


a | lea 


(Figure 12.1.16): 


Figure 12.1.16 


PROOF Both mand M are equal to 1,so 1-A < [fea <1-A. 


D 
COROLLARY 2 


If f(x,y) 20 on D then IJ, fl, y)dd 20. If fly) SO on D then 
{ofl sydd <0. 


To really be sure that the double integral corresponds to the volume, we 
need to know that it is the only operation that has the Addition and Cylinder 
Properties. To make this precise, we introduce the notion of a volume function. 

We suppose f(x, y) is continuous at every point of an open region Dy, and 
consider subregions D of Dy. A volume function for f is a function B which assigns 
a real number B(D) to each subregion D of Dy and has the Addition Property 


B(D) = B(D,) + B(D,) 
and the Cylinder Property 
mA < BID) S MA, 


where m is the minimum and M the maximum value of f on D. 


UNIQUENESS THEOREM 


The double integral \\,, f(x, y) dA is the only volume function for f. That is, 
if Bis a function which has the Addition and Cylinder Properties, then 


B(D) = | [ree y)dA for every D. 
D 
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Given a continuous function f such that f(x, y) = 0 for all (x, y), the function 
V(D) = volume over D 


certainly has the Addition and Cylinder Properties. Thus we are justified in defining 
the volume as the double integral. 


DEFINITION 


Let f(x, y) 2 O for (x, y) in D. Then the volume over D between z = 0 and 
z = f(x, y) is the double integral 


V= [ [ronnaa. 
D 


When f(x, y) is the constant 1, we have 


ae ae 


That is, the area of D is equal to the volume of the cylinder with base D and height 1, 


as in Figure 12.1.17. 
Given any unit of length (say meters), if the height is one meter then the area 
is in square meters and the volume has the same value but in cubic meters. 


Height = 1] - 


Figure 12.1.17 


PROBLEMS FOR SECTION 12.1 


Compute the following double Riemann sums. 


1 2), Gx + 4y) Ax Ay, Ax=4, Ay=4, D:0<x<1, 0<y<1 

2 YY (4 + 2x — Sy)Ax Ay, Ax=4, Ay=4t, D:-25x<52, -I<y<! 
3 55 x? + Ardy, Ax=4, Ay=4, D:-2<x<2, -2<y;<2 

4 YEU + paws Ax=4, Ay=4, D:0<x<2. 0<y7<! 


5 LES Ax ay. Ax=4, Ay=4, D:1l<x<2, 1<y<2 
D 
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6 Dy (cae + siny)Ax Ay. Ax = _ Ay 7 D: 5 <x - O<sy<n 
7 m3 (cosx siny) Ax Ay. Ax = : Ay = o D : Hee Gs yen 
8 YS ve’Ax Ay. Ax =3. Av=l, D:0<5x<52, -2<y <3 
9 5s Pertkeh te: Res Resi De eRBI. Hee peo 
10 sae + 2y)AxAy. Ax=i, Ay={. DiO<x<lb OSy<x 
el See roe Ax=!, Ay=t4, D:0<x<1 xSy<1 
12 ST? +o ARAN: Ax=t. Ay=}, D:-Il<xsl. OSys¥? 
dD 
13 dd ysinxAx Ay. Ax = a Ay = a D:0<x<n. sin’x <y < 2sinx 
14 YS (e* + eX) Ax Ay, Ax=1, Ayv=1, D:-3e5x<3, -xSy<xn 
D 
15 YY 4AxAy, Ax=lo Ayal, Div +r s9 
D 
16 YO -10Axv Ay, Av=t Av=1, D:-3 5x53, ?<y S18 -xX? 
D 
O17 Show that if D is a region with area A and c is constant, then [|, ¢dA = cA. 
O 18 Prove the Constant Rule: 
By ef (x, y) Ax Ay = Owy fin yy Ax Ay, 


¢ [ [res y) dx dy. 


D 


[fe ‘(x y) dx dy 
D 
Oo 19 Prove the Sum Rule: 


ED Sony) + glx yp Ax Ay = YY f(x, y) Ax Ay + VY glx. y) Ax Ay, 
D D D 


[ [res y) + glx y)dx dy = [ [10 y)dx dy + | | a(x, vy) dx dy. 


D D D 


12.2 ITERATED INTEGRALS 


In this section we shall learn how to evaluate double integrals. A double integral 
can be evaluated by two single integrations. The Iterated Integral Theorem gives the 


key formula. 
a bo(x) 
i il I(x, y)d » ts 
a bi(x) 


The iterated integral 
is an integral of an integral of f(x, y). It is evaluated in two stages. First evaluate the 
inside integral 


ba(x) 
g(x) = | Tix y)dy 


bi(x) 
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by ordinary definite integration, treating x as a constant. This gives us a function of x 
alone. Second, evaluate the outside integral 


02, 22 b2(x) 
i a(x) dx -{ | ACA) ar dx 
ay ay b1(x) 


by a second definite integration. 
We shall usually drop the brackets around the inside integral and write the 
iterated integral as 


42 pb2(x) 
[f° se.nayax 


ay ~by(x) 
ITERATED INTEGRAL THEOREM 
Let D be a region 


a,5x< 4a, b,(x) S y S d(x). 


The double integral over D is equal to the iterated integral: 


[fresnaa =f” fs y) dy dx. 


Discussion For a fixed Xo, {224°} f(xo, y) dy is the area of the cross section shown 


in Figure 12.2.1. The Iterated Integral Theorem states that the volume is 
equal to the integral of the areas of the cross sections. 


The proof of the Iterated Integral Theorem is given at the end of this section. 
When using iterated integrals we must be sure that: 


(1) a, S a, and b,(x) S b,(x). 
(2) The differentials dx and dy appear in the right order. 
(3) The outer integral sign has constant limits. 


AWS 


by(xo) 52(xo) y 


Figure 12.2.1 
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While the order of the differentials, dx dy or dy dx, does not matter in a 
double integral, it is important in an iterated integral. The inside integral sign goes 
with the inside differential, and is performed first. 


i (ive r i i 


me 


i 
| 
| 
{ 


: do first : 


do second 


When the region D is a rectangle, there are two possible orders of integration, 
because all the boundaries are constant. Thus there are two different iterated integrals 
over a rectangle. Integrating first with respect to y we have 


a2 pb2 
F(x, y) dy dx, 
ay vb, 
and integrating first with respect to x we have 


be par 
Cee 


by 4a, 


Using the Iterated Integral Theorem twice, we see that both iterated integrals must 
equal the double integral. 


ie f(x, y)dy dx = | [re y)dA 
: : D 

b2 par 
| [ f(x, y) dx dy = J fren dA 


by vay 


Therefore the two iterated integrals are equal to each other. We have proved a 
corollary. 


COROLLARY 


The two iterated integrals over a rectangle are equal: 
a2 bz par 
ii ii S(x,y) dy dx = eh I(x, y) dx dy. 


Discussion Thiscorollary is the simplest form ofa result known as Fubini’s Theorem. 
Remember that by our Permanent Assumption, f(x, y) is continuous on D, 
For an idea of the difficulties that arise when f(x, y) is not assumed to be 
continuous, see Problem 49 at the end of this section. 


There are also other regions besides rectangles over which we can integrate 
in either of two orders, such as Example 5 in this section. 


In the following two examples we evaluate the double integrals which were 
approximated by double Riemann sums in the preceding section. 


12.2 ITERATED INTEGRALS 


EXAMPLE 1 Evaluate [fv dA 
Di 


where D, is the unit square 
O<x<1l, O<yK<l. 


The limits of the outside integral are given by 0 < x < 1, and those of the 
inside integral are given by 0 < y < 1. The iterated integral is thus 


ipl 
[fevaa - | i} xy dy dx. 
x ovo 


The inside integral is 


1 y=l 
{ x*ydy = wy| = 35075 
= 


0 y=0 


1 i 
Then {fy dA = I a tx| = + ~ 0.16667. 
) 
Dy 


i) 


Since D, is a rectangle we may also integrate in the other order, and should 
get the same answer. 


1 i 
[ [rea =f [ xy dx dy. 
f 0 Jo 


1 1 
{ x*ydx =4x3y| =4y. 
0 0 


1 1 

[fey dA = { tydy = »| = + ~ 0.16667. 
0 

Dy 


it) 


The Riemann sums in Section 12.1 were 0.0875, 0.12825. 


EXAMPLE 2 Evaluate |{,,, x*y dA where D, is the region in Figure 12.2.2: 
OC ea4.. sey eae 


The limits on the outside integral are given by 0 < x < 1, and those on the 
inside integral by x? < y < Wes so the iterated integral is 


Figure 12.2.2 
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1 pJx 
[fra =[ f x? dy dx. 
0 vx? 
D2 


fe yale 
i x 2y dy = 5x | = 43 = 4x6, 


= x? 


1 
Ne ydA= [ a (5x? — 4x°) dx = $x* te 
0 


= & ~ 0.05357. 
The Riemann sums in Section 12.1 were 0.05375, 0.04881. 


In many applications the region D is given verbally, and part of the problem 
is to find inequalities which describe D. 


EXAMPLE 3 Let D be the region bounded by the curve xy = 1 and theline y = 3 — x, 
Find fa ders which describe D, and write down an iterated integral 
equal to {f, f(x, y) dA. 


Step 7 Sketch the region D as in Figure 12.2.3. 


Figure 12.2.3 
Step 2. The line and curve intersect where 


x — x)= 1, 


—}3x+1=0, 
(x — $)(x — 2) = 0 
=45, x= 2, 


For 1/2 < x < 2, the curve y = I/x is below the line y = 5/2 — x. Therefore 
D is the region 


Se XS 2; I/x<p<3-x. 


Step 3 The inequalities for x give the limits of the outside integral, and those for y 
give the limits of the inside integral. Thus 


(5/2)-x 
Jf yd =f ie I(x, y) dy dx. 
1/2 
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EXAMPLE 4 Find the volume of the solid bounded by the surfaces z = 0,z = y — x?, 
yell. 


Step 7 Sketch the solid and the region D, as in Figure 12.2.4. 


Figure 12.2.4 


Step 2 Find the inequalities describing the region D. 
This is the hardest step, and gives us the limits of integration. The surfaces 
z = Oand z = y — x? intersect at the curve y = x”. We see from the figure 
that D is the region between the curves y = x* and y = 1, so D is given by 


-lex<1, x sy<1. 


Step 3 Set up the iterated integral and evaluate it. 


1 1 
v=[fy-xda=f{ { y — x? dy dx. 
— {x2 
D 


= (5-1? — x21) — G(x)? — x? - x?) 
1 4 
2 
1 
v=| 4 —x? +4xtdx = 18. 


Multiple integration problems can be solved by a three-step process as shown 
in Examples 3 and 4. 


Step 7 Sketch the problem. 
Step 2 Find the inequalities describing the region D. 
Step 3 Set up the iterated integral and evaluate. 

We can also integrate over a region in the (y, x) plane instead of the (x, y) 
plane. A region D in the (y, x) plane has the form 

b,SySb, ay(y) Sx < ayy), 

as shown in Figure 12.2.5. 

The double integral over D is equal to the iterated integral with dy on the 
outside and dx inside, 
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aly) aly) 


Figure 12.2.5 


[ [te y)dA = le oe I(x, y) dx dy. 


- by Yaiy) 


Some regions, such as rectangles and ellipses, may be regarded as regions in either 
the (x, y) plane or the (y, x) plane (Figure 12.2.6). 


y y 
box) idee eescee 
ay(y) ay) 
| ; 
De eae 
by(x) 
| | 
ay dy ~ x 
(a) Das an (x, y) Region (b) Das a(y, x) Region 
Figure 12.2.6 


EXAMPLE 5 Let D be the region bounded by the curves 
x= yp’, x=yt+2, 
Evaluate the double integral [[, xy dA. 
Step 7 The region D is sketched in Figure 12.2.7. 


Figure 12.2.7 
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Step 2 Find inequalities for D. To do this we must find the points where the curves 
x=y*, x=y+2 
intersect. Solving for y and then x, we see that they intersect at 
(1, —1), (4, 2). 


We see from the figure that D is a region in either the (x, y) plane or the (y, x) 
plane. However, the boundary curves are simpler in the (y, x) plane. D is the 
region 


-lsy<2, y?Sx<y42. 


Step 3 Set up the iterated integral and evaluate. 


2 yr 2 
[ [avaa = { | xy dx dy. 
D pier a 


yt2 


yt2 
| xydx = »| 
y? y2 


= Hy + 2)y — 40") 
ff» dA 
D 


ay? + 2y? + 2y — dy? 
{ ay? + 2y? + 2y — gy° dy = Ye. 
-1 
PROOF OF THE ITERATED INTEGRAL THEOREM For any region D, let B(D) be 
the iterated integral over D. Our plan is to prove that B has the Addition and 
Cylinder Properties, so that by the Uniqueness Theorem B(D) will equal the 
double integral. 


ll 


PROOF OF ADDITION PROPERTY 


Case? Let D be divided into D, and D, as in Figure 12.2.8(a). By the Addition 
Property for single integrals, 


by (x) 


a3 pb2(x) a2 pb2(x) 
- | faydx + | fdy dx 


ay ~b4(x) az ¥ by (x) 


= B(D,) + BID). 


a2 eb2(x) 
B(D) = { f dy dx 


Case 2 Let D be divided into D, and D, as in Figure 12.2.8(b). Then 


a2 pb2(x) 
B(D) = i) { fdy dx 


bi (x) 


a2 pb3(x) b2(x) 
a { | fdy + f aya 


bi (x) b3(x) 


a2 eb3{x) a2 pb2(x) 
- | faydx + | fdydx 


ay ¥ by (x) a; ~b3{x) 


= B(D,) + B(D)). 
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¥ ¥ 
balx) ba{x) 
dD, 
! 
i | ! b | 
1 by(x) | | 
| 
ae 4 f -— 
ay a3 ay 7 ay ay 
(a) (b) 
Figure 12.2.8 


PROOF OF CYLINDER PROPERTY Let m be the minimum value and M_ the 
maximum value of f(x, ') on D For each fixed value of x, 


b2(x) b2(x) 
i mdy < { f(x, y) dy. 


bi(x) bi (x) 


Integrating from a, to a), 


a2 pbo(x) a2 pbo(x) 
| { mdydx < | [ fdy dx = B(D). 
ay Yby(x) ay ¥ by (x) 

a2 


a2 pbo(x) 
But | i m dy dx -| m(b2(x) — by(x)) dx 


bi(x) ay 


a mf (b3(x) — b,(x)) dx = "mA. 


Therefore mA < BD). 
By a similar argument, BID) S MA. 
Since B has both the Addition and Cylinder Properties, 


B(D) = [ron dA. 


D 


The Constant, Sum, and Inequality Rules for double integrals follow easily 
from the corresponding rules for single integrals, using the Iterated Integral Theorem. 


CONSTANT RULE 


[feres y)dA = c| { T(x, yp dA. 
D D 


SUM RULE 


[ [so y) + ex. p)dA = | [re y)dA + { {es yl dA, 
D D D 


12.2 ITERATED INTEGRALS 


INEQUALITY RULE 


If f(x, y) < g(x, y) for all (x, y) in D, 
[[reneas [fac naa 
D D 


PROOF As an illustration we prove the Sum Rule. 


a2 pb2(x) 
[ [r+ saa =f ftedydx 


D a, ¥bi(x) 
a2! b2{x) b2(x) 
-f] fay + | gay ts 
ay bi (x) bi(x) 


az pb2(x) az pb2(x) 
=| fadydx + { i g dy dx 


ay “by (x) bi(x) 
= [ {ra + [feaa. 
D D 


The Iterated Integral Theorem gives another proof that the area of D is equal 
to the double integral of 1 over D. 
By definition of area between two curves, 


A= | (bis) ~ by(x)) dx 
Using iterated integrals, 
a2 pb2(x) 
{fas = i dy dx 
D a; “bi (x) 


a { * (byl) — by(x)) dx = A. 


PROBLEMS FOR SECTION 12.2 


In Problems 1-16, evaluate the double integrals (compare these with the problems from Section 
12.1). 


1 {[ex+anaa, D:0<x<10<y<!l 
D 

2 [f+ 2x-syaa, D:-25x<52,-l<y<1 
D 

3 [for + yaa, D:-2<x<2,-2<y<2 


D 


4 [fas oan, D:0<x<2,0<y<1 


D 


733 


734 


10 


11 


12 


13 


14 


15 


16 
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x 

[ {aa Dilsxs21lsy<2 
y 

D 


[| cos +siny)dA, Di: -aj2Sx<an/205y Sun 


D 


[ftcos siny)dA, Di: -n/2S5xS2n/2,.0S5) Sn 


D 


| [reraa. D:0Sx<2,-2<y;<3 


D 


[ [ermmraa. D:-2<x<2,-2<y)<2 


D 


[foc + anda. D:O0<x<s1d0sy<vx 
D 
[fers+ 38y)dA, D:0SxSskxs<y<l 


D 


{foe + vid, D:-1<x<10<sy<¥x? 
D 

[ [vsinw dA, D:0<Sx<xn,sin?x <y < 2sinx 
D 

[fie + @)dA, Di: -3SxS53,-xS yx 

D 


[aaa Dix? +)? <9 


D 


J { -104a, D:-3<x<3,x7<y<18— x? 


D 


In Problems 17-24, evaluate the iterated integral. Then check your answer, by evaluating in the 
other order. 


17 


19 


21 


23 


1 pl 1 1 
| i (x?y — 3xy? + 5) dy dx 18 [ xy(2y + 1) dy dx 
o*o o*o 
6 8 4 p6 
dy dx 20 3x dy dx 
3¥~2 241 
j2 pni2 1 2 y 
[ [ sin(x + y) dy dx 22 | ) ——s dy dx 
0 *o -1¥9 1 +x 
3 p6 ea 2 al 1 
| i Vx + ydydx 24 ( i dy dx 
a iyo X + J 


In Problems 25-30 evaluate the iterated integral. 


25 


O“sinx 


1 pex rel 
i} i dy dx 26 i i) dy dx 
avo i 
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2 pla—x2 1 pJ1l—x? 
27 { i ydydx 28 { x? + y* dy dx 
iy 0° ~f1—x? 
3 p3y 1 py 2 
29 { { x?y dx dy 30 | | ———— dx dy 
Ov y2 oo ./1 ~— x* 


In Problems 31-38, find inequalities which describe the given region D, and write down an 
iterated integral equal to ff, f(x, y) dA. 


31 The triangle with vertices (0, 0), (5, 0), (0, 5). 

32 The triangle with vertices (1, —2), (1, 4), (5, 0). 

33 The circle of radius 2 with center at the origin. 

34 The bottom half of the circle of radius 1 with center at (2, 3). 

35 The region bounded by the parabola y = 4 — x? and the line y = 3x. 

36 The region above the parabola y = x? and inside the circle x? + y? = |. 

37 The region bounded by the curves x = $and x = I/(1 + y?). 

38 The region bounded by the curves x = 12 + y? and x = y*. 

39 Find the volume of the solid over the region x? + y* < 1 and between the surfaces 


z=0,z=x?. 
40 Find the volume of the solid over the region 
Disgx<2,x<y<x? 
and between the surfaces z = 0, z = y/x. 


41 Find the volume of the solid between the surfaces z = 0, z = 2 + 3x — y, over the 
region0<x<2,0<y<x. 
42 Find the volume of the solid between the surfaces z = 0, z = ,/y — x, over the region 
Osx<i1,xsSy<t. 
43 Find the volume of the solid bounded by the plane z= 0 and the paraboloid 
x? y? 
ol eee nee 
44 Find the volume of the solid bounded by the three coordinate planes and the plane 
ax + by + cz = 1, where a, b, and ¢ are positive. 
45 Show that 


a2 pb2 2 
{ F(x) dy de = (b2 — by) flo) dx. 
46 Show that 
a2 pb2 a2 ba 
i St) + a) dy dx = (b, ~ b,)| Sled + (a, ~ af aly) dy. 


47 Show that 
a2 pb2 a2 bz 
[ [see ay ax = | {s09 as | { 2) ay) 
ay vb a by 
48 Show that 
b palx) 
{ { ydydx = 0. 
ad —g{x) 
_Sy if x 1s rational, 
ai ie Gy) = { 1—y if x is irrational. 
Show that: 
1 1 1 
(a) [ [ fosmavae= [ga =. 
o Jo 0 


(b) Foreach constant yo # 4, the function g(x) = f(x, yo) is everywhere discontinuous, 
so that the iterated integral [3{3 f(x, y) dx dy is undefined. 
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12.3 INFINITE SUM THEOREM AND VOLUME 


The double integral, like the single integral, has a number of applications to geometry 
and physics. The basic theorem which justifies these applications is the Infinite Sum 
Theorem. It shows how to get an integration formula by considering an infinitely 
small element of area. 

An element of area is a rectangle AD whose sides are infinitesimal and parallel 
to the x and y axes. Given an element of area AD, we let 


(x, y) = lower left corner of AD, 
Ax, Ay = dimensions of AD, 
AA = Ax Ay = area of AD. 


AD is illustrated in Figure 12.3.1. 


(x, v) Ax 


(x, ») 


An element of area 
Figure 12.3.1 


INFINITE SUM THEOREM 


Let h(x, y) be continuous on an open region Do and let B be a function which 
assigns a real number B(D) to each region D contained in Dy. Assume that 
(i) B has the Addition Property B(D) = B(D,) + B(D,). 
(ii) B(D) = 0 for every D. 
(iii) For every element of area AD, B(AD) = h(x, y) AA (compared to AA). 


Then B(D) = [ [ms y) dA. 


D 


We shall use the notation 
AB = BAD). 
Given (i) and (ii), the theorem shows that if we always have 
AB = h(x, y) AA (compared to AA) 
then BID) = VY A(x, y) AA. 


D 


The proof is simplest in the case that D is a rectangle. 


12.3 INFINITE SUM THEOREM AND VOLUME 


PROOF WHEN D IS A RECTANGLE Choose positive infinitesimal Ax and Ay and 
partition D into elements of area AD (Figure 12.3.2). Since B has the Addition 
Property, B(D) is the sum of the AB’s. Let c be any positive real number. For 
each AD we have 


AB = h(x, y)AA (compared to AA), 


AB 
na 
AB AB 

Ad SN <4 tS 

AB — cAA < h(x, y)AA < AB + CAA. 


h(x, y), 


Letting A be the area of D and adding up, 
B(D) — cA < VY h(x, y) AA < B(D) + cA. 
D 


Taking standard parts, 


B(D) —cA < { | A(x, y)dA < BUD) + cA, 


D 


so B(D) = [ics y) dA. 


D 


Figure 12.3.2 


The proof in the general case is similar except that some of the elements of 
area AD will overlap the boundary of D and thus be only partly within D. (See Figure 
12.3.3.) The method of proof is to change D to include all instead of part of each AD, 
use hypothesis (ii) to show that the new B(D) is infinitely close to the old one, and 
then show as above that the new B(D) is infinitely close to the double integral 
{[p Ax. y) dA. 

In most applications of the Infinite Sum Theorem, hypotheses (i) and (ii) are 
automatic. To get a formula for B(D) in practice, we take an element of area AD and 
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; o 


Figure 12.3.3 


find an /i(x, y) such that 
AB = h(x, y) AA (compared to AA). 


Our first application is to the volume between two surfaces. 


DEFINITION 


Let f(x, y) S g(x, y) for (x, y) in D and let E be the set of all points in space 
such that 


(x, y) is in D, T(x y) Sz S g(x,y). 
The volume of E is 
v= | fats ~ ¥) — L(x, y) aA. 


V is called the volume over D between the surfaces z = f(x. }) and z = g(x, y) 
(Figure 12.3.4). 


x 


Figure 12.3.4 Volume between two surfaces 
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JUSTIFICATION The part AE of the solid E over an element of area AD is a rectangu- 
lar solid with base AA and height g(x, y) — f(x, y), except that the top and 
bottom surfaces are curved (Figure 12.3.5). Therefore the volume of AE is 


AV & (g(x, y) — f(x, y)) AA (compared to AA). 


By the Infinite Sum Theorem, 


v= [etx — fee ydA. 
D 


Figure 12.3.5 


EXAMPLE 1_ Find the volume of the solid 
O<x<l,- OSy<x, x+y<z<e*” 


Step 71 Dis the triangle shown in Figure 12.3.6. 


Step 2 DistheregionO<x<1,0<y<x. 


Figure 12.3.6 
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Step 3 ae (jer —(x+y)dA 


1px 
= { { etry — (x + y)dydx. 
0v0 


x 
Bs 
| erty = (x +4 y) dy — erty —xXy> we = e2* = phe 3x?, 
c¢) fe) 


1 
Ve [ e?* — eX — 3x7 dx = te? — @. 


EXAMPLE 2. Find the volume of the solid bounded by the four planes 


x = 0, y=0, z=x+y, z=1—-—x-y. 


Step 7 Sketch the planes. We see from Figure 12.3.7 that z= x+y is below 
z=1-x-y}. 


D 
x 
Figure 12.3.7 
Step 2. Find inequalities for the region D. Since the two planes 
z=x+y, z=Il-x-y 

meet at the line 2x + 2y = 1, yot-x, 

D is the region O<x <5, O<sys<t--x. 
Step 3 v= {fa x-—yp)—(x + y)dA= [1 - 2x ~2yaa 

D 


rer 1/2-—x 
= [ | 1 — 2x — 2ydydx. 
0 Jo 
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1j2-x 1j2-x 
[ | =2 = 2ydy=y— dy 7} 
0 0 

3—-x—2x§—x)-G-—xP =4-x4+ x? 


1/2 
es! 1 2 1 
v= | g—X + x°dx = qq. 


EXAMPLE 3 Find the volume of the solid bounded by the plane z = 2y and the 
paraboloid z = 1 — 2x? — y?. 


Step 7 The surfaces and the region D are sketched in Figure 12.3.8. 


=< --- 7 


Figure 12.3.8 


Step 2. The two surfaces intersect on the curve 


2y = 1 — 2x? — y’, 


or solving for y, y= —-1+./2 — 2x’. 


Therefore D is the region 


-1<x<1, —1—./2—2x* << y< -1+./2 — 2x’. 


Step 3 We see from the figure that the plane is the lower surface and the paraboloid 
is the upper surface. 


V= {fa — 2x? ~ y*) — 2ydA 


D 


i -1+/2-2x? 
=| (1 — 2x? — y? — 2y) dy dx. 
-1 


—1-f/2-—2x2 
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-14+ /2—2x? V2-2x? 
(1 — 2x? — y* — 2y)dy -| (2 — 2x? — u*) du 
ei cere er errs 
= 8/2 — x?)3??, 
3 


: : 
2 
V= | ee (Pax Fae 
=1 


Put x = sin0,  ./l — x? = cos@, dx = cos6 d@ (Figure 12.3.9). 


Figure 12.3.9 Jl — x? 
im nj2 
V= “v= [ cos* 6 d0 
3 —nf2 
1 3/1 : 1 m2 
= *y2| cos'dsind + 3[5cos0 sind + 50}) [ 


Ww 
BR) WB 
Nie 

a 

ll 

> 

a 


Answer V= /2n. 


PROBLEMS FOR SECTION 12.3 


Find the volumes of the following solids in Problems 1-8. 


1 O<x<1l, OSyS1l, xypSz<l 

2 O<x<2, OS y<2 wr4+y<z2<8 

3 Osx<l, lsSyps2, xSzsy 

4 O<x<4, OSyS1, xS2z xe’ 

5 Osx<2, OSySx, ySzex 

6 [<xs4, xsps4, ypSxSxy 

7 —-l<x<1l, xSy<S1, xf/y<2sy 

8 O<x<n, ~sinx <y<sinx, —sinx <z<sinx 


In Problems 9-16, find the volume of the solid bounded by the given surfaces. 
9 The planesy =0, x+y=2, z=—-x, z=x 
10 The planesx =0, y=0, 2x4+3y+z2=4 6xt+yp—z2=8 


ll 
13 
15 


17 


18 
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z= x? 4 y?, z=4 12 z=x?4+y? +1, z=2x+2y 

y=0, z=x?+y, z=1 14 x=0, x=y, 2=1l-y 

e+y=9 VP +7=9 16 z=x*4y?, z=2—x7—5P 
2 2 2 


Find the volume of the ellipsoid _ + - + _ a 


Find the volume of the solid bounded by the paraboloid z = x?/a? + y?/b? and the 
plane z = c, where c is positive. 


12.4 APPLICATIONS TO PHYSICS 


In this section we obtain double integrals for mass, center of mass, and moment of 


inertia. 


DEFINITION 


If a plane object fills a region D and has continuous density p(x, y), its mass is 


m= J Joc. naa 


On an element of area AD, the density is infinitely close to p(x, y) (Figure 


12.4.1). Therefore the mass is 


Am ~ p(x, y) AA (compared to AA). 


By the Infinite Sum Theorem m = J, p(x, y) dx dy. 


(Xo, Yo) 


mass An 


Figure 12.4.1 


In Chapter 6 we were able to find the mass of a plane object whose density 


p(x) depends only on x by a single integral, 
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sis [ p(x)(ba(x) — bo) dx. 


Our new formula for mass reduces to the old formula in this case, for by the Iterated 


Integral Theorem, 
m= [ {eos dA 


D 


a2 pb2(x) 
— | p(x) dy dx 


bilx) 


= | p(x)(ba(x) — bi(x)) dx. 


ay 


Now we can find the mass of a plane object whose density p(x, y) depends on both x 
and y instead of on x alone. 


EXAMPLE 1_ Find the mass of an object in the shape of a unit square whose density 
is the sum of the distance from one edge and twice the distance from a second 
perpendicular edge. 


Step 1 The region D is shown in Figure 12.4.2. 


Jy 


Figure 12.4,2 


Step 2 Place the object so the first two edges are on the x and y axes. Then D is the 
region 


Step 3 The density is p(x, y) = y + 2x. 


1 pl 
m =|fy+ 2xdA = [ | y + 2x dydx. 
00 
D 


1 


i 
| y + 2x dy = 4y? + 2x0 = 54 2x. 
0 


0 


1 
m= { 3 + 2xdx = 3. 
0 


DEFINITION 


A plane object which fills a region D and has continuous density p(x, y) has 
moments about the x and y axes given by 


12.4 APPLICATIONS TO PHYSICS 
Mie i | yplx, 9) dA. 
D 


M,= [x00 y) dA. 
D 


M,, and M, are sometimes called first moments to distinguish them from 
moments of inertia (which are called second moments). 


The center of mass of the object is the point (X, y) with coordinates 
{ [xen aa 
M, es 
fl ° 
{{ p(x, y)dA 
D 
ff yp(x, y) dA 
M, D 


mm : 
J {ots aa 


JUSTIFICATION The piece of the object on an element of area AD has mass 


x= 


y= 


Am = p(x, y) AA (compared to AA). 
A point mass m at (x, y) has moments 
M,. = ym, M, = xm. 
Therefore the piece of the object at AD has moments 
AM, x yAm & yp(x, y) AA (compared to AA), 
AM, = x Am & xp(x, y) AA (compared to AA). 


The double integrals for M, and M, now follow from the Infinite Sum 
Theorem. 


An object will balance on a pin at its center of mass (Figure 12.4.3). The 
center of mass is useful in finding the work done against gravity when moving the 


Center of mass 


Figure 12.4.3 
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object. The work is the same as if the mass were all concentrated at the center of mass, 
and is given by 
W = mgs 


where s is the distance the center of mass is raised and g is constant. 


EXAMPLE 2 A triangular plate bounded by the lines x = 0, x = y, y= 1 has 
density p(x, y) = x + y. Find the moments and center of mass. 


Step 7 Sketch the region D, as in Figure 12.4.4, 


Figure 12.4.4 


Step 2. We see from the figure that D is the region 


O<x<1, xSySt. 


Step 3 Set up and evaluate the iterated integrals for the mass m and moments M,, 


and M,. 
1 pl 
m= I} +ydA =| { x + ydydx. 
Ovx 
D 


1 
| x+ ydy=x+4— 3x’. 
x 


1 
m= [ x +4 — 3x? dx =}, 
0 


1 al 
M,= [{yx+yda={[ [ yx + y? dy dx. 
Ovx 
D 


=a ee a Pr 
M, =| 9X43 - 4x7 dx = 3%. 


4 
0 
1 pl 
: { [ + y)dA = | i x? + xydydx 
0 vx 
D 


= 
Il 


12.4 APPLICATIONS TO PHYSICS 


The answers are M, = a M,= 
es M, 5/24 5 
m 1/2 12 
_ M, 9/24 9 
a 1p 12 


The point (x, ¥) is shown in Figure 12.4.5. 


EXAMPLE 3 The object in Example 2 is lying horizontally on the ground. Find the 
work required to stand the object up with the hypotenuse of the triangle on 
the ground (Figure 12.4.6). 


We use the formula W = megs. 


From Example 2, m = 4. We must find s. 


ae : 5 9 ; 
Ss = minimum distance from lal to the line x = y. 
bie 5 \2 g \2 
s = minimum value of z = Je _ 5 + (s = a . 
28 106 
= 2s? aun 
Ory Ea ia 
dz 7 28 cans 
dx 12} 2 
dz 28° 7 
fee ee ag 


The second moment, or moment of inertia, of a point mass m about the origin 
is the mass times the square of the distance to the origin, 


Figure 12.4.5 Figure 12.4.6 
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I = m(x* + y?). 


The moment of inertia is related to the kinetic energy of rotation. A mass m moving at 
speed v has kinetic energy 


KE =3 mv’. 


Hence if m is rotating about the origin with angular velocity w radians per second, 


its speed ist = w./x? + y? and 
KE = 4n(m/x? + y’)? = Ho?. 


Thus moment of inertia is the rotational analogue of mass. 


DEFINITION 


Given a plane object on the region D with continuous density p(x, y), the 
moment of inertia about the origin is 


I= ie yx? + y*) dA. 
D 


JUSTIFICATION On an element of volume AD, the moment of inertia is 
Al = (x* + y?)Am & p(x, yx? + y?)AA (compared to AA). 
The integral for J follows by the Infinite Sum Theorem. 


EXAMPLE 4 Find the moment of inertia about the origin of an object with constant 
density p = 1 which covers the square shown in Figure 12.4.7: 


1 a 1 1 1 
7SxXS3, —zSySsz. 


1/2 1/2 
{Jo + y*)dA = | x? + y? dy dx. 
-—1/2 
D 


— 1/2 


— 
ll 


1/2 1/2 
i x? + ypdy=x*y+ | =x 4h. 


— 172 


x? + aedx =}. 


by 
l| 
t 
es, = 
: N 
nN 


Figure 12.4.7 
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PROBLEMS FOR SECTION 12.4 


In Problems 1-10, find (a) the mass, (b) the center of mass, (c) the moment of inertia about the 
origin, of the given plane object. 


1 —asxSa —bsy<b, plx,y)=k 


2 O<x<a, O<y<b, p(xy)=k 

3 O<x<1, x <y<1, piyy=k 

4 O<x<a, OS y<bx, plxy)=k 

5 O<x<2, x<y<2x, p(xy)=xt+ytl 
6 O<x<sl, OSy<x, p(x,y)=x—-y 

7 O<xsl1, Os y<>x’, p(x») = /x + /y 
8 1sxS2, xSySx?, p(x,y) = 1//xy 

9 OS<x<52, e*syse, p(xy=1 

10 -I<x<1, O<ysI//1+x’, p(x yay 

11 Find the mass of an object in the shape of a unit square whose density is the sum of the 
four distances from the sides. 

12 Find the mass of an object in the shape of a unit square whose density is the product of 
the distances from the four sides. 

13 An object on the triangle 0 < x < 1,0 < y < x has density equal to the distance from 
the hypotenuse y = x. Find the amount of work required to stand the object up (a) on 
one of the short sides, (b) on the hypotenuse. 

14 An object in the shape of a unit square has density equal to the distance to the nearest 
side. Find the mass and the amount of work needed to stand the object up on a side. 

15 An object on the plane region —1 <x < 1, x? <y <1 has density p(x, y)=1+x+ /y. 
Find the mass and the work needed to stand the object up on the flat side. 

16 An object on the unit square 0 < x < 1,0 < y <1 has density p(x, y) = ax + by +c. 
Find the mass and center of mass. 

oO 17 The moment of an object of density p(x, y) in the region D about the vertical line x = a 
is defined as 
Myaaa= { [ (x ~ adele y)dA. 
Show that D 
My, .=a= M,—a-m 
where M, is the moment about the y-axis and m is the mass. 
oO 18 The moment of inertia of an object in the region D of density p(x, y) about the point 


P(a, b) is defined as 


Ip= { [oes ne = a)? + y= by dA. 
Show that D 
Ip = 1 — 2aM, — 2bM, + m(a? + b?) 


where J is the moment of inertia about the origin, M, and M, are the first moments, and 
m is the mass. 


12.5 DOUBLE INTEGRALS IN POLAR COORDINATES 


A point with polar coordinates (8, r) has rectangular coordinates 


(x, y) = (r cos 6, * sin 6). 
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DEFINITION 
A polar region is a region D in the (x,y) plane given by polar coordinate 
inequalities 
as0< Bp, a(0) <r < b(0), 


where a(@) and b(@) are continuous. To avoid overlaps, we also require that 
for all (0,1) in D, 


O<O0<2n and Or. 


The last requirement means that the limits « and f are between 0 and 2z, 
while the limits a(@) and b(@) are =O. Figure 12.5.1 shows a polar region. 
The simplest polar regions are the polar rectangles 


2<0< 8, asr<b. 


We see in Figure 12.5.2 that the 0 boundaries are radii and the r boundaries are 
circular arcs. 


olar region 
Figure 12.5.1 SS DOME ASE! 


A polar rectangle 
Figure 12.5.2 2 e 


The polar rectangle 
a2<0<68, Os<r<b 


is a sector of a circle of radius b (Figure 12.5.3(a)). 
The polar rectangle 


O<0< 2a, Ox<r<b 


is a whole circle of radius b (Figure 12.5.3(b)). 


Less trivial examples of polar regions are the circle with diameter from (0, 0) to 
(0, 5), 
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b 
b 
B 
a 
(a) sector: (b) circle: 

asé<B,0<r<b 0<@<27,0<r<sb 

Figure 12.5.3 
O0O<@<z, O<r<bsin8g, 

and the cardioid 0<6 <2z, 0<r<1-4cos@. 


Both of these regions are shown in Figure 12.5.4. 

We shall use the Infinite Sum Theorem to get a formula for the double 
integral over a polar region. In the proof we take for AD an infinitely small polar 
rectangle. 


POLAR INTEGRATION FORMULA 
Let D be the polar region 
as6<f, a(8) <r < b(6). 
The double integral of f(x, y) over D is 
B pb(0) 
{| I(x, y)dA = [ F(x, yyr dr dé 
D 


a(@) 


B pb(@) 
=| i f(r cos 8, r sin )r dr dO. 


a(@) 


Notice that in the iterated integral for a polar region we do not integrate 
f(x, y) but the product of f(x, y) and r. Intuitively, the extra r comes from the fact 
that a polar element of area is almost a rectangle of area r A@ Ar (see Figure 12.5.6(b)). 


PROOF We shall work with the rectangular (@,1r) plane. Let C be the region in the 
(9, r) plane given by the inequalities 
ax<O0<f, a(é) <r < b(8). 


Thus C has the same inequalities as D but they refer to the (8, r) plane instead 
of the (x, y) plane. D and C are shown in Figure 12.5.5. 


| J f(x, y) dx dy = J | f(x, yr dO dr. 


D G 


We must prove that 
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(a) circle: (b) cardioid: 
OS 6S 7,05 rs bsiné 0< 65 27,05 r<1+cos é 


Figure 12.5.4 


a: 


Figure 12.5.5 


Our plan is to use the Infinite Sum Theorem in the (0,1) plane. Assume 
first that f(x, ) > 0 for all (x, y) in D. 


For any (0,r) region C, corresponding to a polar region D, in the (x, y) 
plane, let 


B(C,) = [ [406 y) dx dy. 
Dy 


Then B has the Addition Property and is always =O. Consider an element of 
area AC in the (0,r) plane with area A@ Ar. AC corresponds to a polar 
rectangle AD in the (x, y) plane. As we can see from Figure 12.5.6, AD is 
almost a rectangle with sides r AO and Ar and area r A@ Ar. 


The volume over AD is almost a rectangular solid with base of area r AQ Ar 
and height 


f(x, ») = f(r cos 6, r sin 0). 
Therefore B(AC) = f(x, yr AO Ar (compared to A@ Ar). 
By the Infinite Sum Theorem 


B(C) = [ [7 weded 
Cc 


and by definition B(C) = [10 y) dx dy. 


D 
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(r, 6) 


(a) 


(b) Exaggerated scale 
Figure 12.5.6 


Finally we consider the case where f(x, y) is not always positive. Pick a 
real constant k > 0 such that f(x, y) + k is always positive for (x, y) in D. 
By the above proof, 


[{ vow + Waxdy = [{ Ye») + radar, 


D c 
{{ kdxdy = {] kr dO dr. 
D Cc 


When we use the Sum Rule and subtract the second equation from the first, 
we get 


[ [70 y)dx dy = {{ S (x, y)r dé dr. 


D Cc 


In a double integration problem where the region D is a circle or a sector 
of a circle, it is usually best to take the center as the origin and represent D as a polar 
rectangle: 


EXAMPLE 1 Find the volume over the unit circle x? + y* < 1 between the surfaces 
z= Oandz = x’. 


Step 1 Sketch D and the solid, as in Figure 12.5.7. 
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Figure 12.5.7 


Step 2. Dis the polar region 0 < @ < 22,0 <r 


<t. 


2n 1 2n 1 
Step3 V= {] x? dx dy = i x?r dr dQ = { { (r cos 6)*r dr dO 
F o Jo 0 Jo 
2n 
-[ (ic. 1 cos? @ dr dé 


2a 
-[- 4cos? 6. d@ = 44 sin 0 cos 8 + 10] = 7/4. 
0 _10 


For comparison let us also work this problem in rectangular coordinates 
We can see that it is easier to use polar coordinates 


Dis the region -1 <x <1, -/l—-x* sys fl — x? 
J1—x2 1 
ale *dxdy = fo [ x? dy dx = [ 2x7. /1 — x? dx. 
Ji= x2 - 


We make the trigonometric substitution shown in Figure 12.5.8 


x=sing, /1—x*?=cos¢, 


dx = cos ddd. 


Figure 12.5.8 V1 — x? 


Then @ = —n/2 atx = —land ¢@ =27/2 at x = 1,s0 


nj2 
v= | 2 sin? d cos? 6 dh 


—nj2 


ate (' “7 (| teste dp 


2 
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= ie: 3(1 — cos? 26) db 


—n{2 
n 


= i 4(1 — cos? u) du 


t 


= 7(u — 4 cosusinu — 5u = 7/4. 
Yu - 3 4 


EXAMPLE 2 Find the mass and center of mass of a flat plate in the shape of a semi- 


circle of radius one whose density is equal to the distance from the center 
of the circle. 


Step 1 The region D is sketched in Figure 12.5.9. 


Step 2. Take the origin at the center of the circle and the x-axis as the base of the 
semicircle. D is the polar region 0< 0<2,0<r< 1. 


Step 3 The density is 


A(x, y)= f/x? +y? =r. 
mz pl 
m= |{ x? + ydA -| { rer dr d@ 
o Jo 
D 
inn 2 dr dO [5 = 
—— r ‘g ——4 _ ==, 
o Jo o 3 3 
nr opl 
Mo | rors =| ip rsin@-r-rdr dé 
= [ r? sin 0 dr dO -| zsin 6 dé = 3. 
My, a1) x4 y?dA= [Ff reosa-r-rdr ao 


=f { 79 cos O dr d8 = | 4cos 0d0 = 0. 
0 40 0 


R M, _ M 3 
Answer m= 3? X= res 0, y 5) 77 


The point (x, ) is shown in Figure 12.5.10. 


Figure 12.5.9 Figure 12.5.10 
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EXAMPLE 3 Find the moment of inertia of a circle of radius b and constant density 
p about the center of the circle. 


Step 1 Draw the region D (Figure 12.5.1 {). 


Figure 12.5.11 


Step 2 Put the origin at the center, so D is the polar region 


0<@S8 2a, O<rsb. 


Xe) 
2n pb 
I= [f e-0 + y*)dA= [ [ pr? +r dr d 
0 vo 
D 


2n pb 2n 
=p| [ Pardo =p ib* do 
o vo 0 


7 pb*n 
 <e 


Step3 x*4+y? =P, 


PROBLEMS FOR SECTION 12.5 


In Problems !~16, find the volume using polar coordinates. 
1 xr+y<l, O<z<6 


2 e+yre<sl, O<zsx74 5 

3 voy<4 0<z<x4+2 

4 e+ ye s4, O<2<s fe +¥ 

5 ty <9, x+y? <2 <9 

6 x+y S25 O<2se"” 

7 bsx?ty? <4, OP 4 ppl sz s(x? + yy? 
8 Pe 4 ys9> te ea Ses 

9 S421, CH ye Joe, OSS a5 
10 —$2y2) OSyey/e= et ve se p79 
5) n4s0< 27, OSr<1, O<z2<Pr 

12 O<O0<n/6, 1Srs2, OSz25 9-7 

13 O<O0<n, OSr<2sind, OSz2<r 

14 O0<O0<n/2, OSr<cost, P<z<Pr? 

15 Os0<2a, OS rsd, 0528 77034 2r0 
16 0S0s2n, 0USrse. 0=2< 2/7 
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17 Find the volume of the solid over the cardioid r = 1 + cos @ between the plane z = 0 
and the cone z =r. 

18 Find the volume of the solid over the cardioid r = 1 + cos @ between the paraboloids 
z=r and z= 8 — 7. 

19 Find the volume of the solid over the circle r = sin 0 between the plane z = 0 and the 
hemisphere z = ./1 — 7°. 

20 Find the volume of the solid over the circle r = 2 cos 6 between the plane z = 0 and the 
cone z= 2 —7. 

21 Find the volume of the solid over the polar rectangle « < @ < B,a <r < b, between the 
plane z = 0 and the cone z = r. 

22 Find the volume of the portion of the hemisphere 0 < z < ./1 — 7? over the polar 
rectangle a < 0 < B,a <r <b (assuming b < 1). 

23 A circular object of radius b has density equal to the distance from the outside of the 
circle. Find (a) the mass, (b) the moment of inertia about the origin. 

24 A circular object of radius b has density equal to the cube of the distance from the center. 
Find (a) the mass, (b) the moment of inertia about the origin. 

25 Find the moment of inertia about the origin of a circular ring a <r < b,0 <0 < 2h, 
of constant density k. 

26 Find the moment of inertia of a circular object of radius b and constant density k about 


a point on its circumference. (The center can be put at (0, b), so the object is on the polar 
region0 <r < 2bsin0,0<0<7z,) 


27 An object has constant density k on the circular sectorO<x<1,0<y<./1 — x? 
Find (a) the center of mass, (b) the moment of inertia about the origin. 

28 An object of constant density k covers the cardioid r < 1 + cos 6,0 < 6 < 2z. Find 
(a) the center of mass, (b) the moment of inertia about the origin. 

29 An object of constant density k covers the region inside the circle r = 2bsin @ and 
outside the circle r = b. Find (a) the center of mass, (b) the moment of inertia about the 
origin. 

30 An object of constant density k covers the polar region 


0<0< 72/2, O<r< bsin 20. 


Find (a) the center of mass, (b) the moment of inertia about the origin. 


O 31 (a) Use polar coordinates to evaluate ; { 5 gene? dy dx. 


([ ea} 


wn 


(b) Show that | [- cays = [f e-* dx 


= 20 


(c) Now evaluate the single integral e-* dx, 


12.6 TRIPLE INTEGRALS 
A closed region in space, or solid region, is a set E of points given by inequalities 
a, Sx Sap, b,(x) < y S b,(x), c(x%, y) Sz S c,(x, y) 
where the functions b,(x), b,(x) and ¢,(x, y), c2(x, y) are continuous. 


The boundary of E is the part of E on the following surfaces: 


The planes x = a,;, X = 4@y. 
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The cylinders y = b,(x), y = (x). 

The surfaces z = ¢,(x, y), 2 = C(x, y). 

The simplest type of closed region is a rectangular solid, or rectangular box, 
ad; Sx Sap, b, Sy S bo, Cy S ZS Cp. 


Figure 12.6.1 shows a solid region and a rectangular box. 


ay 


x ae 


Region in space Rectangular box 
Figure 12.6.1 


An open region in space is defined in a similar way but with strict inequalities. 
As in the two-dimensional case, the word region alone will mean closed region. 


PERMANENT ASSUMPTION 


Whenever we refer to a function f(x, y, 2) and a solid region E, we assume that 
T(x, y, 2) is continuous on some open region containing E. 


The triple integral If T(x, y, 2) dx dy dz 
E 


is analogous to the double integral. 
The first step in defining the triple integral is to form the circumscribed 
rectangular box of E (Figure 12.6.2). This is the rectangular box 


a, Sx Sa), B, Sy Bz, C,Sz85C), 
where B, = minimum value of b,(x), 
B, = maximum value of 53(x), 


C, = minimum value of ¢,(x, y), 
C, = maximum value of ¢,(x, y). 
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Figure 12.6.2 The circumscribed rectangular box 


Our next step is to define the triple Riemann sum. Given positive real 
numbers Ax, Ay, and Az, we partition the circumscribed rectangular box of E into 
rectangular boxes with sides Ax, Ay, and Az (Figure 12.6.3). The partition points of 
this three-dimensional partition have the form 

(Xz Vis Zm)s O<k<a, O<!l<p, O<megq. 
The triple Riemann sum of f(x, y, z) Ax Ay Az over E is defined as the sum 


LLLSO y,z)AxAyAz= YYY  f(%y, Ys Zm) Ax Ay Az. 


(Xk, ¥1,Zm) in E 


When we replace Ax, Ay, Az by positive infinitesimals dx, dy, dz we obtain an infinite 
triple Riemann sum 


YY f(x, », 2) dx dy dz. 
E 


A partition 
Figure 12.6.3 x 
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LEMMA 
For all positive infinitesimals dx, dy, and dz, the triple Riemann sum 


VDD S06 yz) dx dy dz 
E 


is a finite hyperreal number and therefore has a standard part. 


We are now ready to define the triple integral (see Figure 12.6.4), 


Figure 12.6.4 x 


DEFINITION 


Given positive infinitesimal dx, dy, and dz, the triple integral of a continuous 
junction f(x, y, z) over E is 


{ff f(x,y, 2) dx dy dz = st}Y’YY f(x, y, z) dx dy dz}. 
E 
B 


We shall now briefly state some basic theorems on triple integrals, which 
are exactly like the corresponding theorems for double integrals. 


INDEPENDENCE OF ax, dy, AND az 


The value of ie T(x, y, 2) dx dy dz does not depend on dx, dy, or dz. 


We shall usually use the notation dV = dxdydz for the volume of an 
infinitesimal dx by dy by dz rectangular box, and write 


TJ fO.y,z)dV for ne I(x, y, 2) dx dy dz. 


ADDITION PROPERTY 


If E is divided into two regions E, and E, which meet only on a common 
boundary then 


J See) aE HI T(x. yz) dV + I] f(x, y, 2 dV. 
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ITERATED INTEGRAL THEOREM 
If E is the region 


ap SxS ay; b,(x) S y S b,(x), cy(x, y) S Zz S c,(x, y), 
b2(x) pert, ” 
then [fa x, y,z)dV= ie { I I (x, y, z) dz dy dx. 
bi(x) Yei(x,y) 


If the region E is a rectangular box 
a, Sx <a, b, Sy <by, ¢,Sz5c, 


there are six different iterated integrals over E, corresponding to six different orders 
of integration. Here they are (in “alphabetical” order). 


a2 pbr per a2 pc2 pbo 
(1) | | flx.y,2z)dzdydx (2) { (So ayaeds 
a, ~by 4¥ey a, “Cy by 


ba par C2 b> C2 ar 
i fl,y,2)dzdxdy (4) i) fl, y,2)dx de dy 
b by Cy 


ai 


c2 par pbo c2 pbz par 
i ) T(x y, 2) dy dx dz (6) { I | I(x, y, 2) dx dy dz. 
C1 ay vb; C1 by vay 


The Iterated Integral Theorem shows that each of these iterated integrals is equal to 


the triple integral 
Ff sosnaar. 
E 


EXAMPLE 1 Evaluate {ff xy?z? dV where E is the rectangular box 
E 
0<x <2, Osy<l, O<z<4. 


There are six iterated integrals which all have the same value. We compute 
one of them, and then another to check our answer. 


2 1 4 
FIRST SOLUTION [ [fore ar = | i { xyz? dz dy dx. 
d 0 v0 ¥0 
E 


The inside integral is 


The second integral is 


The final answer is 
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Ae AM 
SECOND SOLUTION {fy xy?z3 av= | [ xy?z? dy dx dz. 
0 ¥0 40 


E 


The inside integral is 
1 


1 
| xy?ziidy = bye = ixz?, 
0 19) 


The second integral is 
2 


2 
| ixz3 dx = je] = $73. 


ra) i) 


The final answer is 


4 4 
4,3 1,4 256 _ 128 
[tet de = aet| = = 


0 0 


Triple integrals can be evaluated by iterated integrals. 


EXAMPLE 2 Evaluate fy y + zdV where E is the region shown in Figure 12.6.5, 
E 


O<x < 7/2, Os ySXsinx, OS<zS ycosx. 


Figure 12.6.5 
n{2 psinx pycosx 
SOLUTION {ff ytzdV= } | | y + zdz dy dx. 
0 0 0 
E 
We first evaluate the inside integral. 
y¥ COS x ycos x 
| ytzdz=yz+t+ HF = y? cos x + 3? cos? x. 
0 0 


Now we evaluate the second integral. 


sin x 


sinx 
{ y? cos x + 4y? cos? x dy = 4y3(cos x + $cos*x)| 
0 


0 
= $sin* x(cos x + 4 cos? x). 


Finally we evaluate the outside integral. 
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0 


nf2 
[fo +zdV= { 4 sin? x(cos x + 4.cos? x) dx 
E 
ni 
= | 4(1 — cos? x)(cos x + 4.cos? x) sin x dx 
oO 


0 
= | —4(1 — u’)(u + 3’) du 
i 


1 
1 1,2 3. 1,4 19 
=] 3(u + 3u u 2u") du = Tg0- 
0 


COROLLARY 


The volume of a region E in space is equal to the triple integral of the constant | 
over E as illustrated in Figure 12.6.6, 


v= ff 


Figure 12.6.6 


PROOF E is the solid over the plane region D given by 
a, xq, b,(x) S y S ba(x) 


between the surfaces z = c,(x, y) and z = c,(x, y). By definition of the 
volume between two surfaces, 


v= [fobs — exts, yaa. 
D 
Using the Iterated Integral Theorem, 
a2 pbx) pe2(x,y) 
{{fav=| [ dz dy dx 
a, ¥bi(x) “c1(x,y) 
E 
a2 pb2(x) 
=[{ Ceey) - exes sayde =v. 
a, ¥bi(x) 


We now come to the Infinite Sum Theorem for triple integrals, which is, 
again, the key result for applications. 
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We shall use Ax, Ay, and Az for positive infinitesimals. By an element of 
volume we mean a rectangular box AE with sides Ax, Ay, and Az (Figure 12.6.7). 
The volume of AE is 


AV = Ax Ay Az. 


Figure 12.6.7 An element of volume 


INFINITE SUM THEOREM 


Let h(x, y, z) be continuous on an open region E, and let B be a function which 
assigns a real number B(E) to each region E contained in Ey. Assume that: 


(i) Bhas the Addition Property. 
(ii) B(E) = 0 for every E. 
(ili) For every element of volume AE, 


B(AE) = h(x, y, z) AV (compared to AV), 
Then B(E) = If h(x, y, z) dV. 
E 


Here are some applications of the triple Infinite Sum Theorem. Perhaps 
the simplest physical interpretation of the triple integral is mass as the triple integral 
of density. 


DEFINITION 


The mass of an object filling a solid region E with continuous density p(x, y, Z) is 


m= {I p(x, y, z) dV. 
E 


JUSTIFICATION At every point of an element of volume AE the density is infinitely 
close to p(x, y. Z), so the element of mass is 


Am = p(x, y, Zz) AV (compared to AV), 
(See Figure 12.6.8.) By the Infinite Sum Theorem, 


n= VJ p(x, y, z) AV. 
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Figure 12.6.8 


EXAMPLE 3 Find the mass of an object in the unit cube 
0<=%=1, OS ysl, OS 7<1 


with density p(x, y,z7)=x+y4+z. 


m= {{fery + sav 


E 
1 1 1 

-| {fj x + y+ z2dzdydx 
0 Jo Yo 


i I 1 
- | { xy tidydx = | eS Sux Ss. 
0 40 


8] 


An object in space has a moment about each coordinate plane. 


DEFINITION 


If an object in space fills a region E and has continuous density p(x, y, z), its 
moments about the coordinate planes are 


M,, = SVJ zp(x, y, z) dV. 
My. = J] frat naav 


M,, = J | fave y, z) dV. 


The center of mass of the object is the point (x, y, Z), where mis mass and 


M M M,., 
y= —, ye = 7=—. 
m m m 


JUSTIFICATION A point mass m has moment M,, = mz about the (x, y) plane 
(Figure 12.6.9). In an element of volume AE, the object has moment 


AM,, = z Am = zp(x, y,z) AV (compared to AV). 
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Figure 12.6.9 


By the Infinite Sum Theorem, 
My = If zplx, y, Zz) dV. 
E 


EXAMPLE 4 An object has constant density and the shape of a tetrahedron with 
vertices at the four points 


(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), 
Find the center of mass. 


Step 1 The region is sketched in Figure 12.6.10. 


Figure 12.6.10 


Step 2 The region E ts the solid bounded by the coordinate planes and the plane 
X + y+ 2= 1 which passes through (1, 0,0), (0, 1.0), (0,0, 1). Solving for 


z, the plane is 
z=!1-—-N-y. 
This plane meets the plane z = 0 at the line 1 -x —y =O, ory =1—-x. 
Therefore E is the region 
O<x<l, O<ysIt-~y, O<zs!1-x-y. 
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Step 3 Let the density be p = 1. 


1 1-x pi-x-y 
m={{[av= | | i dz dy dx 
0 Yo 0 
E 
1 pl-x 
- | | 1—x — ydydx 
o Yo 


M,.= [[{xav 


th Mie: 1/24 1 
ve n~(«i/6 4 
Similar] j= i z= ! 
The center of mass is 
a 111 
(x, ¥, Z) = aa} 


An object in space has a moment of inertia about each coordinate axis. 


Intuitively, the moment of inertia about an axis is the analogue of mass for rotations 
about the axis. 


DEFINITION 


If an object in space fills a region E and has continuous density p(x, y, 2), 
its moments of inertia about the coordinate axes are 


I, = {ff (y? + 2 )a(x, ys z) dV, 
E 

I= {ff (x* + 27)p(x, y, 2) dV, 
E 

i= {J (x? + y?)o(x, y, z) dV. 
E 


JUSTIFICATION A point mass m has a moment of inertia about the x-axis of 


l= (y? + 27)m. 
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On an element of volume AE, the object has moment of inertia 
Al, = (y? + 27) Am = (y? + 2*)p(x, y,2) AV (compared to AV). 
The triple integral for I, follows by the Infinite Sum Theorem. 


EXAMPLE 5 Find the moments of inertia about the three axes of an object with 
constant density 1 filling the cube shown in Figure 12.6.11, 


O<x<a, Osy<a, O<z<a. 


n= [[Jore P42dV= inne + 2? dzdy dx 
=| [ ay” + 3a° dy dx = q* dx = 40°, 
0 


1, 


Similarly, 


Figure 12.6.11 


PROBLEMS FOR SECTION 12.6 


In Problems !-8, evaluate the iterated integral. 


1 pd 3 
1 {| [ svedzay dy 
0 42 1 
2 p2 pd 
2 | { } (x — 2p + 42) dzdydx 
i 0 = 
4 1 a 
3 [ [, ii (By? + 627) dz dy dx 
2 p2 
‘ 
0 
1 
2 ak 
0 
1 
ot 
0 
7 an cs " /ya dz dy dx 
0 40 
2 


tf cOosx pysinx 
8 ( { [ (x + 2z)dzdydx 
0 


0 0 


fre e**? dzdydx 
i x? + yz) dz dy dx 


2x?z dz dy dx 
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In Problems 9-16, evaluate the triple integral. 


9 [ [e+ 2nav. E:0Sx52,1<y<3,2<72<4 
E 

10 J {fever ar, E:0<x<10<y<10<z<1 
E 

1 [[ fa + sear, E:0<x<10,0<y<x7,0<2<xy 
E 

1 2 3 

12 {ff Sat dv, E:l<x<el<y<xl<z<x 
pi lx voz 

13 LL etertar, E:-1<x<ix<y<1,x<z<y 
E 

14 [ff xersar, E:1<x<2,0<y<Inx,0O<z<y 
E 

15 [[[ Ver or ear, E:0<x<10<y<x,yp<z<2y 
E 

16 {{fav. E:0Sx<Lx?<y<x,x*yszsx/yp 
E 


In Problems 17-26, find (a) the mass, (b) the center of mass, (c) the moments of inertia about the 
three coordinate axes, of an object with density p(x, y, z) filling the region E. 


7 E:0<5x<1,0<y<10<5z<1 P(x, y,Z) = x + 2p + 3z 

18 E:0<x<10<y<10<z<1  p(x,y2=x?+y42? 

19 E:0sx<lO<sysx,0<z<x4+y p(x, y, Z) = 2 

20 E:-1Sx<i,x?<y<1x?<z<1 ) p(x,y,z)=z 

21 E:0<xS1,0Sy<1fxysz<s1 p(x, y, Z) = xyz 

22 E:0sx<1lxsysix<z<y p(x, y, 2) = 10 

23 E is the tetrahedron with vertices at (0, 0, 0), (a, 0, 0), (0, b, 0), (0, 0, ¢), p(x, y, 2) = k. 

24 Eis the tetrahedron with vertices (0, 0, 0), (1, 0, 0), (1, 1, 0), (1,0, 1), p(x, y,z) =x + y4z. 
25 E is the rectangular box0 <x <a@,0< y<b,0<2z<e, p(x,y,z)=k. 

26 E is the rectangular box -a<x<a,-bsy<b,-cSz<oplx,y,2J=k. 


CYLINDRICAL AND SPHERICAL COORDINATES 


In evaluating triple integrals it is sometimes easier to use cylindrical or spherical 
coordinates instead of rectangular coordinates. 
A point (x, y, z) has cylindrical coordinates (0, r, z) if 


x =rcos@, y =rsin 0, 222: 


That is, as we see in Figure 12.7.1, (6,1) is a polar coordinate representation of (x, y), 


and z is the height above the (x, y) plane. 
The name cylindrical coordinates is used because the graph of the cylindrical 
coordinate equation r = constant is a circular cylinder as shown in Figure 12.7.2. 
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r = constant 


Cylindrical coordinates x 
Figure 12.7.1 Figure 12.7.2 
DEFINITION 


A cylindrical region is a region E in (x,)',Z) space given by cylindrical co- 
ordinate inequalities 


“2<0< BP, ado) < r < b(0), c(0,r) Sz S c,(0,r), 


where all the functions are continuous. To avoid overlaps we also require 
that for (0,6, 2) in E, 


O<@<2n and O<r. 


A cylindrical region is shown in Figure 12.7.3. 
The simplest kind of cylindrical region is the cylindrical box 


x4S08 8B, axrx<b, Cy S250). 


This is a cylinder whose base is a polar rectangle and whose upper and lower faces 
are horizontal, as in Figure 12.7.4. 


c(8, r) 


A cylindrical region A cylindrical box 
Figure 12.7.3 Figure 12.7.4 
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The cylinder box 
0<0<2z, O<r<b, ¢, S78 ey 


is a cylinder whose base is a circle of radius b and whose top and bottom faces are 
horizontal (Figure 12.7.5). 
The cylindrical box 


0<6@ < 2z, a<r<b, 6p 2 ZZG 


OS@S2r 
po oe O<r<b 
OS zS cy 


Peete tate 


— 


Figure 12.7.5 Figure 12.7.6 


To get a formula for the triple integral over a cylindrical region E, we use 
the Infinite Sum Theorem but take for AE an infinitely small cylindrical box instead 
of rectangular box. 


CYLINDRICAL INTEGRATION FORMULA 
Let E be the cylindrical region 
asO0sf, a(0) < r < b(8), c,(0,r) < z < c,(0,r). 
The triple integral of f(x, y, z) over E is 


b(8) c2(8,. if, 
JJ Jers y,z)dV = [ | { S (x, y, zr dz dr dO. 
a(0) Jc1(6,r) 


To evaluate the triple integral we substitute 
f(x, y, Z) = f(r cos 6, r sin 6, z). 
This is like the Polar Integration Formula but has an extra variable z. 
In the iterated integral we do not integrate f(x, y, z) but the product of f(x, y, z) and r. 
PROOF Let C be the region in the rectangular (0, r, z) space given by 
aesO< Ff, a(@) <r < bi), c;(0,r) Sz < c,(8,7). 
The region C is shown in Figure 12.7.7. 
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v4 


Figure 12.7.7 


We must prove that 


J J J Laer a de IJ fx, y, zr dO dr dz. 


Assume first that f(x, y, z) > 0on E. For any (0, r, z) region C, corresponding 
to a cylindrical region E,, define 


B(C,) = [[f16 y, 2) dx dy dz. 
E, 


B has the Addition Property and is 20. An element of volume AC in the 
(8, 1", z) space has volume AO Ar Az. As we can see from Figure 12.7.8, AC cor- 
responds to a cylindrical box AE. AE is almost a rectangular box with sides 
r AO, Ar, and Az, and volume r A@ Ar Az. 


At any point of AE, f has value infinitely close to 
S (x, y, 2) = f(r cos 8, r sin 0, 2). 
Therefore B(AC) = f(x, y, z)r A@ Ar Az (compared to A@ Ar Az). 


Figure 12.7.8 
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By the Infinite Sum Theorem 
B(C) = {ff T(x, y, z)r dé dr dr, 
c 
and by definition 


BIC) = [fr y, z) dx dy dz. 
E 


The general case where f(x, y, z) is not always positive is dealt with as in the 
Polar Integration Formula proof. 


When integrating over a solid region E whose base is a circle or polar 
rectangle, it is often easier to use cylindrical instead of rectangular coordinates. 


EXAMPLE 1. Find the moment of inertia of a cylinder of height h, base a circle of 
radius b, and constant density 1, about its axis. 


Step 7 Draw the region as in Figure 12.7.9. 


Step 2. The problem is greatly simplified by a wise choice of coordinate axes. Let 
the z-axis be the axis of the cylinder and put the origin at the center of the 
base. Then the region £ in rectangular coordinates is 


“b= x <b, SP? aps JP x7; Of rsh 


and in cylindrical coordinates is 
O<@0<2n, O<rsb, OSz<h. 


Step 3 The problem fooks easier in cylindrical coordinates. 


x? + y? = 4? 


2n pb ph 
L, =| {fo + Pyav={ { rr dz dr dO 
_ 0 0 40 


2n pb ph 
- | { { r° dz dr dé 
0 Jo vo 
nb*h 


2 pb 5 an] F 
= har db = { —b*hd@ = ——. 
a ae 


Figure 12.7.9 
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EXAMPLE 2 Find the center of mass of a cone of constant density with height h 


and base a circle of radius b. 
The region is sketched in Figure 12.7.10. 
Put the origin at the center of the base and let the z-axis be the axis of the 
cone. £ is the cylindrical region 
O0O<@ 2n, O<r<b, O<2<h—F, 


Step 7 
Step 2 


Step 3 Let the density be 1. 


a 2n pb ph—hrjb 
Ne *| dz dr do 
0 Jovo 
2n 
-[ [. fhe | dr 
2m] nb*h 


L:. 5 1 
= b3 0 = 2 = 
ie | b*h i 4 rau h do 3 


a 3b 


e2n ab ph—hrib 
= {ff zaVa | { | zr dz dr dO 
7 0 0 40 


2n pb a\ 2 
= [ I, a _ “) dr dO 
1 2n 
=o [ [ Pos ra dO 


2 3 4 
= he [5 erate 
0 3 


2 4b? 
i mb? h? 

=i? | <p? d= 
x i. pe 5 


Since the cone is symmetric about the z-axis, ¥ = 


Mt os, ot 1 
Z7=—“=-h, X,¥,Z) = (0,0,-h}. 
- m 4" (72) | | 


The point (X, ?, Z) is shown in Figure 12.7.11. 


(%, ¥, 2) 


1 
\ 
| 
i 
1 
1 
A ht 
ie | 


fan 


Figure 12.7.10 Figure 12.7.11 


12.7 CYLINDRICAL AND SPHERICAL COORDINATES 775 
To express a point P(x, y, z) in spherical coordinates we \et p (rho) be the 
distance from the origin to P, let @ be the same angle as in cylindrical coordinates, 


and let @ be the angle between the positive z-axis and the line OP. Note that @ can 
always be chosen between 0 and z. 


(¢, 9, p) 


Figure 12.7.12 


We see from Figure 12.7.12 that a point (x, y, z) has spherical coordinates (0, @, p) if 
x = psin ¢ cos 6, y=psin ¢sin 8, z= pcos ®¢. 


The graph of the equation p = constant is a sphere with center at the origin (hence 
the name spherical coordinates). The graph of ¢ = constant is a vertical cone with 


vertex at the origin. The graph of @ = constant is a half-plane through the z-axis. 
These surfaces are shown in Figure 12.7.13. 


p = constant 


@ = constant 


* 6 =constant 
Figure 12.7.13 


DEFINITION 


A spherical region E is a region in (x, y, 2) space given by spherical coordinate 
inequalities 
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x Ss 0 Ss Zs B (8) s o = B(8), c,(8, p) s p Ss c,(0, $)s 
where all the functions are continuous. To avoid overlaps we also require 
that for (0, @, p) in E, 

O0< 0 2n, O0<¢d<a, OS p. 


A spherical box is a spherical region of the simple form 
a, SAS 4%, Bi SOS Ba, Cy Spey. 


The 0-boundaries are planes, the @-boundaries are portions of cone surfaces, and 
the p-boundaries are portions of spherical surfaces. Figure 12.7.14 shows a spherical 


box. 


Figure 12.7.14 A spherical box 


The spherical box 
0<6@S 2n, O0<o<x, O<psc 
is a sphere of radius c with center at the origin. 
The spherical box 
0<6@< 22, O<¢d<f, OSpse 
is a cone whose vertex is at the origin and whose top is spherical instead of fat. 
(See Figure 12.7.15.) 


Cone with spherical top 


Figure 12.7.15 


12.7. CYLINDRICAL AND SPHERICAL COORDINATES 777 


Another important example is the spherical region 
0<@<2z, 0<¢<n2/2, 0<p<ccos#¢, 


which is a sphere of radius $c whose center is on the z-axis at z = 3c (Figure 12.7.16). 


Figure 12.7.16 


When integrating over a solid region E made up of spheres or cones, it is 
often easiest to use spherical coordinates. 


SPHERICAL INTEGRATION FORMULA 


Let E be a spherical region 
a S054, B<PS BO), c,(6.6) <p <c,(0, 4). 
The triple integral of f(x, y, z) over E is 


x2 pB2(8) pc2(d,6) 
| f i fl, y,2)4V = | Sc, y,2)p? sin } do dd dd. 
E a 


Bx(9) “ c1(,8) 


In practice we make the substitution 
f(x, y, Z) = f(p sin ¢ cos 0, p sin @ sin 8, p cos d) 


before integrating. 


PROOF Let C be the region in the rectangular (0, d, p) space which has the same 
inequalities as E. We prove 


[[[ fooserdxavae = [{ f s0s.y, 200? sin 6 a9 46 dp, 
E Cc 
As usual we let f(x, y, z) > 0 on E and put 


B(C,) = {ff T(x, y, z) dx dy dz. 
Ey 


Consider an element of volume AC. As we see from Figure 12.7.17, AC 
corresponds to a spherical box AE. AE is almost a rectangular box with 


sides 
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p sin @ Ad 
Figure 12.7.17 Spherical Element of Volume 


Ap, pAd, psind Ad 
and volume p? sin ¢ AO Ad Ap. 
Thus B(AC) = f(x, y, z)p? sin @ A@ Ad Ap (compared to A@ Ad Ap). 


By the Infinite Sum Theorem 
B(C) = {}f S(x,y, 2)p? sin b dO d¢ dp, 
c 


and by definition 


B(C) = al ST (x, p, 2) dx dy dz. 
E 


The triple integral for volume, 


|) a 


gives us iterated integral formulas for volume in rectangular, cylindrical, and 
spherical coordinates. 


ar pbx) pers.) 
Rectangular V= } i dz dy dx. 
ay 


by(x) Yer(xy) 


B ab(6) pc2(8, r) 
Cylindrical v= | i ene One 
a a(@) c,(6,r) 
a2 pB8) pe2(b,d) 
Spherical = [ p* sin b dp dd do. 
vx, YB (0) ¥c (6,9) 


The rectangular formula is really equivalent to the double integral for the 
volume between two surfaces. Similarly, the cylindrical formula is equivalent to the 
double integral in polar coordinates for the volume between two surfaces. 


12.7 CYLINDRICAL AND SPHERICAL COORDINATES 


On the other hand, the volume formula in spherical coordinates is something 
new which is useful for finding volumes of spherical regions. 


EXAMPLE 3 Find the volume of the region above the cone ¢ = f and inside the 
sphere p = c. 
The region, shown in Figure 12.7.18, is given by 


0<¢@8 22, 0<¢SfB, O<pe. 
2n pB pe 
v=| [ [ e?sin bap do ao 
0 Jo Jo 
2n B 3 
= { Ssin odo ao 
o Jo 


ra 3 
-[ (1 cos f)S a0 = (1 cos B)c?. 


1) 


Figure 12.7.18 


EXAMPLE 4 A sphere of diameter a passes through the center of a sphere of radius 
b, and a > b. Find the volume of the region inside the sphere of diameter a 
and outside the sphere of radius b. 


Step 7 The region is sketched in Figure 12.7.19. 


p=acos¢ 


Figure 12.7.19 
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Step 2. We let the z-axis be the line through the two centers and put the origin at 


the center of the sphere of radius b. The two spheres have the spherical 
equations 


p =acos ¢, p=b. 
; b 
They intersect at cos¢@ =-. 
a 
Thus E is the region 


b 
0<@<2z, 0 < @ S arccos -, b<p<acos@. 
a 


arccos (b/a) acos@ 
Step 3 y=fo I I p’ sin & dp df d0 


arccos (b/a) { 
a [ 3 @ cos? ¢ — b?)sin ¢ d¢ dd. 


Put u = cos ¢,du = —sin ¢ dd. Then 


2n pbja 1 
= | — =(a3u? — b) dud 
o 4 3 


1 2n J ae 4 
=3], { avue — b° du dé 
= bja 


ae 3 bt b* 
- Be aed pads Pe eee . 
=3 4 + re dé “(« 4b° + 3 =| 


EXAMPLE 5 Find the mass of a sphere of radius c whose density is equal to the 
distance from the surface. The sphere is shown in Figure 12.7.20. 


Figure 12.7.20 


Put the center at the origin. The sphere is then given by 
O< 6S 2n, O< <x, Ope. 
The density at (0, @, p) is 
density = c — p. 
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The mass is 


m= is i if (c — p)p” sin b dp do dé 


0 040 
2n pt 
_ C3 _ ale 
=| ( Ge je*) sin 6d ao 
2 
ria I, 4 
= _~c* sin ¢ dg dé 
[ [ qsetsin ao 
2n{ é 
= =c* d@ = —c* 
e 


PROBLEMS FOR SECTION 12.7 


In Problems 1-6, evaluate the integral using cylindrical coordinates. 


1 


[fv + y?zdV, Eis the cylinder x? + y? <1,0<z<2 
E 


{{f +zdV, Eis the cylinder x? + y?<9,0<z<6 
E 


[[Je+re E is the cone x? + y?<1,0<z2<1—/x? + y? 
[[f4+ Vea, E is the cone x? + y* <1,./x* + y?<z<1 
E 


fy (x + y)zdV, Eistheregion0<x<2,0<y</4—-x7,0<z2<x?+y 
E 


{ {| ew E is the region 1 < x* + y? <4,0 <z < Jx| 

E 
Find the mass of an object in the shape of a cylinder of radius b and height h whose 
density is equal to the distance from the axis. 
Find the mass of an object in the shape of a cylinder of radius b and height h whose 
density is equal to the distance from the base. 


Find the mass of an object in the shape of a cone of radius b and height h whose density 
is equal to the square of the distance from the axis. 


Find the mass of an object in the shape of a cone of radius b and height h whose density 


is equal to the sum of the distance from the base and the distance from the axis. 


Find the center of mass of an object of constant density filling the region above the 
paraboloid z = x? + y* and below the plane z = 1. 
Find the center of mass of an object of constant density filling the region 

x? + y? <b, O<z</x? + y?. 
Find the moment of inertia of an object of constant density k in the cylinder 0 < r < b, 
—c < zc, about the x-axis. 


Find the moment of inertia of an object of constant density k in the cylindrical shell 
as<r<b,—-c<z<c, about the z-axis. 


781 


782 


15 


16 
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Find the moment of inertia of an object of constant density k in a cone of radius b and 
height 4 about its axis. 

Find the moment of inertia of an object of constant density & in a cone of radius b and 
height 4 about a line through its apex and perpendicular to its axis. 


In Problems 17--24, evaluate the integral using spherical coordinates. 


17 


18 


19 


20 


21 


22 


23 


{ff x? + y? +2? dV, Eis the sphere x? + 5? + 2? <b? 
£ 

I] / fx + y? + 22dV, Eis the sphere x? + y? + 2? < b? 
Je 


ies dV, Eis the sphere x? + y? + 2? <1 


lA 


{fl z? dV, Eis the sphere x* + y? + 2? 
E 

ff [ zav. E is the sphere p < 2bcos ¢ 
Je 

iy (x? + y? + 2?)9? dV,  E is the intersection of the spheres p < 2bcos ¢,p <b 
E 


{ff zx? + y? + 2 dV, E is the region above the cone ¢ = « and inside the sphere 
E 
p=b 


{fl SaaF dV, Eis the spherical shell a < p < b 
E 


Find the volume of the spherical shell a < p < b. 

Find the volume of the spherical box «, S$ 0S %.,8; <@S B,¢, Sp Se. 

Find the volume of the region above the cone ¢ = f and inside the sphere p = bcos @. 
Find the volume of the spherical region 0 < 6 < 27,05 @< 2,05 p sing. 


Find the mass of an object in the shape of a sphere of radius c whose density is equal 
to the distance from the center. 


Find the mass of a spherical shell «@ < p < b whose density is equal to the reciprocal 
of the distance from the center. 


Find the moment of inertia of a spherical object of radius b and constant density k 
about a diameter of the sphere. 


Find the moment of inertia of a spherical shell a < p S b of constant density k about 
any diameter. 


A hole of radius @ is bored through a sphere of radius b, and the surface of the hole 
passes through the center of the sphere, a = 3b. Find the volume removed. 


A hole of radius a is bored through a cone of height 4 and base of radius b, and the 
axis of the cone is on the surface of the hole (a < 4b). Find the volume removed. 


Find the center of mass of a hemisphere of constant density and radius b. 
Find the moment of inertia of an object of constant density k in the ellipsoid 


x y Zz 
Stigtaye=! 
a bt 3 c 


about the z-axis. Hint: Change variables x, = x/a, y; = y/b, z, = 2/c and use spherical 
coordinates. 


EXTRA PROBLEMS FOR CHAPTER 12 


EXTRA PROBLEMS FOR CHAPTER 12 


1 


12 


13 


14 


15 


16 


17 


21 


22 


Compute the Riemann sum 
YY x? — yp? Ax Ay, Ax 4,Ay £, D: 2 S.x.5 2) =2:5 yp S.2. 
D 
Compute the Riemann sum 


DD x? — \/y Ax Ay, Ax =4,Ay=4, D:-1<x<10<y<1-x2 
D 


Evaluate { { x? —ydA, D:-2<x<2,-2<y<2. 
D 


Evaluate { [ x? - Jy dA, D:-1<x<1,0<y<1-x?. 


D 
2/2 p/i— 2x2 

Evaluate { { xdy dx, 
0 


— J/1-2x2 
1 1 
Evaluate [ { ye* dy dx, 
0 J /x 


Find the volume of the solid over the region —1 < x < 1,0 < y < 1 — x* and between 
the surfaces z= 0,z = 1 — y. 
Find the volume of the solid over the region x? + y? = 4 and between the surfaces 
z=Oandz=y?+x+42, 
Find the volume of the solid 1 < x < 2,0 <y <inx, y/x <z < I/x. 
Find the volume of the solid x? + y? < 1,x7j3 <z<1. 
Find the volume of the solid bounded by the planes 
xX SK, Bea Z=x+y, z=x4+2. 

Find the volume of the solid bounded by the cylinders 

xe+yal, x?4¢2=1. 
Find the mass, center of mass, and moment of inertia about the origin of the plane 


object : 
O<x<n, O<y<sinx, plxyy=k. 


Find the mass, center of mass, and moment of inertia about the origin of the plane object 
O<x<1, x<ySl, plxyy)=x?y. 
A circular disc filling the region x? + y? < 1? has density p(x, y) = y?. Find the mass, 
center of mass, and moment of inertia about the origin. 
A semicircular object on the region 
-reix<yr, O<y< Jy? _ x2 
has density p(x, ¥) = y. Find the work required to stand the object up on its flat side. 
Using polar coordinates, find the volume of the solid 
e+yr<9 yS<z<x45. 
Find the volume of the solid over the region 0 < + < 3 + cos@ between the plane 
z = Oand the cone z = r. 
Find the volume of the solid over the circle 0 < r < a between the plane z = 0 and the 
surface z = I/r. 
Find the mass and the moment of inertia about the origin of a semicircular object 
0<r<1,0 < 6 <7z whose density is p(r, 0) = ré. 
A plane object covers the circle 0 < r < a and its density depends only on the distance 
r from the center, p(r, @) = f(r). Show that the center of mass is at the origin. 


nm pmf2 pl 
Evaluate i { i zsinx + zcos ydzdydx. 
QO v0 0 
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23 


25 


26 


27 


29 


30 


31 


32 


33 


34 


35 
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x2 


1 5 
Evaluate [ ( 


0-0 


Evaluate the triple integral 


2 


y 
[ x+y4+zdzdydx. 


0 


[{[>e@ay, eisxs4isysxisesy, 
E 


An object has constant density & in the region 
E:0<xslOsys1l1-x,0<5z xy. 

Find its center of mass and its moments of inertia about the coordinate axes. 
Use cylindrical coordinates to evaluate ie z dV, where E is the region inside the cylinder 
x? + y* = 1 which is above the plane z = 0 and within the sphere x* + y? + 27 = 9. 
An object of constant density k has the shape of a parabolic bowl 

O<O0<2n, OSrsbh rPSzsPrte. 
Find its center of mass and its moment of inertia about the z-axis. 


Use spherical coordinates to evaluate the integral 


iff xtytzdy, 
E 
E is the spherical octant 


ert+yt2<l, OSx, OSy, Oz 
A spherical shell a < p S b has density equal to the distance from the center. Find its 
mass and its moment of inertia about a diameter. 
Prove that the double Riemann sum ya T(x, y) dx dy is finite whenever /(x, y) is 
continuous, D is a closed region, and dx, dy are positive infinitesimals, 


Suppose a plane object is symmetric about the x-axis, that is, it covers a region D of 
the form 


Diasx Sb, —g(x) Sy S g(x) 
and has density p(x, y) = p(x, — y). Prove that the center of mass is on the x-axis. 


The moment of inertia about the x-axis of a point in the plane of mass m is J, = my~ 
Use the Infinite Sum Theorem to show that the moment of inertia about the x-axis 
of a plane object with density p(x, y) in the region D is 1, = |p p(x, y)y? dA. 


The kinetic energy of a point of mass m moving at speed v is KE = 4mv?. A rigid object 
of density p(x, y) in the plane region D is rotating about the origin with angular velocity 
@ (so a point at distance d from the origin has speed wa). Use the Infinite Sum Theorem 
to show that the kinetic energy of the object is 


KE= [| 4q*(x? + y*)o(x, y) dA = 4007. 
‘D 

Suppose a plane object is symmetric about the origin; that is, it fills a polar region 

OsrSg(0), -—x <@ <n, such that g(@ + m) = g(8), and its density has the property 

plr, 8) = p(r, 8 + 2). Show that the center of mass is at the origin. 

Use the Infinite Sum Theorem to show that if D is a polar region of the form a < r < b, 

a(r) < 0 < f(r), then 


{f ff y)dA = [ ee: Ow dé dr. 


D a ¥2x(r) 


